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Abstract 

We report on the search for a charged partner of the X(3872) in the decay B — > X^K, X^ — > 
J/i/jir^TT^, using 213 million BB events collected at the T(45) resonance with the BABAR detector 
at the PEP-II e^e" asymmetric-energy storage ring. The resulting product branching fraction 
upper limits are B[B^/B^ X^K^, X^ J/ipTr^ir^) < 5.8 x lO^^ and S(B=^ ^ ^"^^5: 
X^ J/tpir^TT^) < 11 X 10^® at the 90% confidence level. All results are preliminary. 
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1 INTRODUCTION 



Since the discovery of the X(3872) by the Belle Collaboration PQ, there have been experimental 
confirmations from the CDF [2], DO and BABAR 111 Collaborations. Numerous theoretical 
explanations have been proposed for this high-mass, narrow- width state decaying into J/ipTT^Tr~ . 
The possibilities include a charmonium state [HI, a meson molecular state jTj, and a hybrid 
charmonium state |H]- The Cornell potential model fHI predicts a charmonium state, previously 
unseen, with quantum numbers n'^^^^Lj=l^D2 and J^*-" = 2 . Also this charmonium state 
should be a narrow width state with a 3.830 GeV mass and have a large radiative transition rate, 
X(3872) 7Xci) which has not been observed by Belle jj. Since the measured X(3872) mass is 
very close to the mass threshold, another attractive possibility is a molecular model which 

is a bound state of mesons. If such states exist, then bound states of charged and neutral mesons 
or charged molecular states are plausible. 

The TT~^TT~ mass distributions from the X(3872) decay measured by Belle and BABAR 
both peak near the kinematic upper limit and may be consistent with the decay of it^-k~. 
However, due to limited statistics a spin-parity analysis has not been performed. If indeed, the 
observed decay is X(3872) J/ipp^, then a charged partner, X(3872)='= — > J/t/^p^, may exist. 
Assuming the X charged partner is a member of an isotriplet and isospin is conserved in the 
B decays, the decay rate of B X^K should be twice that of i? ^ X^K. This would make 
experimental detection of the X^ quite favorable. To test this conjecture, the BABAR collaboration 
has performed a search, presented in this paper, for the B-meson decays, B^/B*^ X^K^ and 
B^ X^K^, where X^ J/V'7r=^7r°. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis correspond to a total integrated luminosity of 193 fb~^ taken 
on the T(45) resonance, producing a sample of 213.2ib2.3 million BB events (N^^). Data were 
collected at the PEP-II asymmetric-energy e~^e~ storage ring with the BABAR detector, which is 
described in detail elsewhere TO^. The BABAR detector includes a silicon vertex tracker (SVT) 
and a drift chamber (DCH) in a 1.5 Tesla solenoidal magnetic field to detect charged particles 
and measure their momenta and energy loss (dE/dx). Photons, electrons, and neutral hadrons 
are detected in a CsI(Tl)-crystal electromagnetic calorimeter (EMC). An internally reflecting ring- 
imaging Cherenkov detector (DIRC) provides particle identification information that is complemen- 
tary to that from dE/dx. Penetrating muons and neutral hadrons are identified by resistive-plate 
chambers in the steel flux return (IFR). Preliminary track-selection criteria in this analysis follow 
previous BABAR analyses |11| and a detailed explanation of particle identification (PID) is given 
elsewhere |11| . jl2j . 

3 ANALYSIS METHOD 

This analysis commences with charged and neutral track selections. Each charged track candi- 
date is required to have at least 12 DCH hits and transverse momentum greater than 100 MeV/c. 
If it is not associated with a Kg decay that track candidate must originate near the nominal beam 
spot. 

A charged-kaon or -pion candidate is selected on the basis of dE/dx information from the SVT 
and DCH and the Cherenkov angle measured by the DIRC. An electron candidate is required to 
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have a good match between the expected and measured energy loss (dE/dx) in the DCH, and 
between the expected and measured Cherenkov angle in the DIRC. The ratio of EMC shower 
energy to DCH momentum, and the number of EMC crystals associated with the track candidate 
must be appropriate for an electron. A muon is selected on the basis of energy deposited in the 
EMC, the number and distribution of hits in the IFR, the match between the IFR hits and the 
extrapolation of the DCH track into the IFR, and the depth of penetration of the track into the 
IFR. 

A photon candidate is identified from energy deposited in contiguous EMC crystals summed 
together to form a cluster which has total energy greater than 30 MeV and a shower shape consistent 
with that expected for an electromagnetic shower. 

The intermediate states in the neutral B^/B^ — > J/i/jTr^ir^K^, and charged B^ J/tpTT^ir^Kg 
decay modes used in this analysis are J/if) ^ e~^e~, J/if) ^ /^^/^~; 77) and Kg vr+vr". 

They are selected to be within the mass intervals 2.95 < M(e^e^) < 3.14, 3.06 < M < 3.14, 

0.119 < M(77) < 0.151, and 0.4917 < M(7r+7r") < 0.5037 GeV/c^. The e+e" mass interval is 
larger than that for /i^/x^ in order to recover events in which part of the energy was carried away by 
bremsstrahlung photons. The orientation of the displacement vector between the decay vertex 
and the J/i/j vertex in the lab is required to be consistent with the Kg momentum direction. 

The search for B signal events utilizes two kinematic variables [TJ: the energy difference A^^ 
between the energy of the B candidate and the beam energy in the T(4S) rest frame; and the 

beam-energy-substituted mass tties = \/ (^b)^ ~ (Pb)^' where is the reconstructed momentum 
of the B candidate in the T(45') frame. Signal events should have m-Es ~ "^B) where rriB is the 
nominal mass of the i3-meson, and \AE\ ~ 0. 

Before the data were analyzed, the selection criteria were optimized and fixed separately 
for the charged and neutral modes using a Monte Carlo (MC) simulation of signal and known 
backgrounds. The number of reconstructed MC signal events (nf^^) and the number of recon- 
structed MC background events (n^^) in the signal-box were used to estimate the sensitivity ratio 
n™'^/(a/2 + Y^n™*^) jl3j . where a, the number of standard deviations of significance desired, was 
set to 3. Note the maximum of this ratio is independent of the unknown signal branching fraction. 
This ratio was maximized by varying the selection criteria on AE, rriEs-, the X{J/il)Tr^Tr^) mass, 
the Kg{TT~^TT^) mass, the Kg decay length signficance, the 77 mass, and the particle identification 
for electrons, muons and charged kaons. When there are more than one candidate (on average 
there were 1.3 candidates/event) per event, the candidate with the smallest absolute AE value was 
chosen. All the following plots are displayed with one candidate per event. 

The AE and tues data distributions, after applying the optimized cuts for the neutral and 
charged B modes, are shown in Figs. ^ and El respectively. A clear signal peak is observed at zero 
in the AE distribution and near 5.279 GeV/c^ in the mES distribution. The other feature in the 
AE plots is a wide peak near 0.2 GeV/c^ which is due to i? ^ J/tpK* events combined with a 
random pion track. The rectangular area (signal-box region) bounded by |mES ~ ^b| < 5 MeV/c^ 
and \AE\ < 20 MeV was found to be optimal to select signal events. Choosing events in the 
signal-box region and applying a mass cut of 0.67 < M(7r±7r°) < 0.78 GeV/c^ to select the mass 
region, the K^tt^tt^ mass distributions are shown in Fig. |21for the charged and neutral B modes. 
There are clear signal peaks for ivr['(1270) K^p^ and K^{\21f)) —> Kgp^ corresponding to the 
decays, B^ — > J/ipKf^ and B^ — >■ J/ipKi, previously observed by Belle ^1]. In Fig. El the dashed 
histogram background estimates are obtained using the rriEs sideband region, 5.24 < rriEs < 5.26 
GeV/c^. The number of observed Ki events are consistent with the Belle measurements. 

The J/t/jTT^TT^ mass spectra from the neutral and charged B modes are shown in Fig. [l] No 
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charged decay signal, J/t/jir^TT^, is evident at 3.872 GeV/c^. The mass spectra have back- 

grounds that peak near 3.7 GeV/c^ and have a step near 4.0 GeV/c^. Prom MC studies we find 
the peak near 3.7 GeV/c^ is due to V'(3686) — > J/iJjt^'k decays where one pion is exchanged with a 
random vr*^. The step near 4.0 GeV/c^ is caused by i? — > J /if^Ki, Ki pK decays. 



4 RESULTS AND SYSTEMATIC UNCERTAINTIES 



To extract an upper hmit for — > J/ipir^Tr^, requires a search for a signal in the J/ipTr^ir^ 
mass, mEs, and AE distributions. A signal from B — > X^K, X^ — > J/ipTT^ir^ should pro- 
duce signal peaks in all three distributions. The peaking background from non-resonant, B — >■ 
J/iIj-jt^tt^K, would produce peaks in the ruEs and AE distributions and a flat J/ipir^-iT^ mass 
distribution near 3.872 GeV/c^. The combinatoric background will not create peaks in any of the 
three distributions and should produce a mES distribution whose shape can be parametrized by 
an ARGUS function .17 . To estimate the number of signal events (n^), we count the number 
of observed events {tLoi^s) in the signal region and subtract the estimated number of combinatoric 
background events (ricomb) and the estimated number of peaking background events (upeak)- 

The number of observed events, riobs, is obtained by counting the number of events satisfying, 
\mES - msl < 5 MeV/c^, \AE\ < 20 MeV/c^, and M (J/^vr^vrO) - 3.872MeV/c^ < 12 MeV/c^. 

The number of combinatoric background events, ricomb, is extracted from the thes distribution 
obtained after requiring |A£;| < 20 MeV/c^, and M (J/V'vr±7r°) - 3.872MeV/c^ < 12 MeV/c^. 

The J/i/jTT^TT^ signal band has a 24 MeV/c^ wide mass window. These niEs distributions for the 
neutral and charged B modes are separately fit with the sum of a signal Gaussian function and 
an ARGUS function. The histograms with the fits are shown in Figs. Inland [7| for the neutral and 
charged B modes, respectively. The resulting ARGUS function is integrated over the niES range, 
\mES — ^b\ < 5 MeV/c^, to produce ricomb- The error acomb is obtained from the fit error on the 
normalization of the ARGUS function. The resulting values for ricomb and CTcomb are listed in Table 
1. 

The number of peaking background events, ripcak, is extracted from the rriEs distribution ob- 
tained after requiring |A^| < 20 MeV/c^, and 48 < M (J/V'7r=t7r°) - 3.872MeV/c^ < 72 MeV/c^. 

This J/ipTT^TT^ sideband has a 48 MeV/c^ wide mass window and is twice the mass range of the 
signal band. These ruEs distributions for the neutral and charged B modes are separately fit with 
the sum of a signal Gaussian function and an ARGUS function. The ruEs histograms with the 
fits are shown in Figs. El and [3 for the neutral and charged B modes, respectively. The estimated 
number of peaking background events , ripcak, is calculated by counting the number of events in 
the \mES ~ rriBl < 5 MeV/c^ ruEs region, subtracting the number of combinatoric events obtained 
from integrating the ARGUS function over the same range, \mEs — rnsl < ^ MeV/c^, and finally 
dividing the result by two. Note the Gaussian has a width that was fixed to a value that determined 
from a fit to the rriEs distribution obtained using both the J /ipir vrO si gnal band and the J /ipir n 
sideband. The error apeak is obtained by adding in quadrature the Poisson errors on the number of 
events in \mES ~ < 5 MeV/c^ and the fit errors on the normalization of the ARGUS function. 
The resulting values for rip^ak and a peak are listed in Table 1. 

The total background (nf,) is the sum of the peaking and combinatoric backgrounds and its 
error ((7^) combines in quadrature the errors from the peaking and combinatoric backgrounds. The 
backgrounds and their errors are summarized in Tabled 
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Table 1: Efficiencies, number of signal-box events, and estimated number of background events 
(peaking, combinatoric, total) for the neutral and charged B decays. 



Mode 


e 




'^peak i ^peak 




rib ± <^b 




10.65% 


87 


31.2 ± 8.0 


70.6 ± 6.3 


101.8 ± 10.2 




8.50% 


31 


0.6 ±4.7 


27.0 ±4.0 


27.6 ± 6.2 



The efficiencies (e) for the processes, B^/B^ X^K^, J/i/jtt^tt^ and B^ X^K^, 

X^ J/ipir^T:^ are determined by MC simulation using an X^ signal with zero width, mass 
3.872 GeV/c^ and a model consisting of the sequential isotropic two body decays B X^K, 
X^ — > J/tpp^ and -k^t:^. Efficiencies are corrected for the small differences between data and 

MC by using well-understood control samples where results from data and MC are available. These 
corrections are applied to PID, neutral detection, and tracking efficiencies. The final efficiencies for 
each mode are listed in Tabled 

The systematic errors include uncertainties in the number of BB events in the data sample, 
the secondary branching fractions, the MC statistics, the decay model for the generated events, 
the background parametrization, the particle identification, the charged particle tracking, and 
the vr'^ reconstruction. The individual uncertainties are given as percentages in Table |21 The 
secondary branching fractions include B{J/ip — > e^e~ , p^ =0.1181 =b 0.0014 and B{Kg — 
7r"'"7r~)=0. 686 lb 0.0027. The decay model uncertainty is estimated by comparing the efficiencies for 
phase space and different decay models with J^*^ = 1++ and J^*^ = 2 . 

Table 2: Percentage Systematic Errors from the neutral and charged B decay modes. 



Systematic Errors(%) 






No. of BBevents 


1.1 


1.1 


Branching fractions 


5.3 


5.3 


MC statistics 


2.1 


2.3 


MC decay model 


1.1 


3.0 


Bkgd sideband width 


0.8 


1.9 


Particle ID 


5.0 


5.0 


Tracking tt* 


1.4 


1.4 


Tracking K"^ 


1.4 




Tracking ifg 




2.6 


Tracking J/tp — > e^e~,/i^^~ 


1.8 


1.8 


tt" correction 


3.2 


3.2 


TOTAL (asys) 


8.8 


9.7 



The background parametrization uncertainty is estimated by varying the background sideband 
width, refitting the mss distributions, and recalculating the number of events. The uncertainties in 
particle identification, charged tracking efficiency and vr'^ reconstruction efficiency are determined 
by studying control samples JJ. The total fractional errors {(Jsys) listed at the bottom of Table [21 
are determined by adding the individual contributions in quadrature. 

The probability distribution of the signal events is modeled as a Gaussian with a mean [ng) and 

sigma (cTs). For each i3-decay mode the mean is Ug = Uobs - nb and the sigma is ag = J nobs + o'l 
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X 1 + a'^yg . The systematic error is added in quadrature and scales the errors on Uobs and Ub by 
the same fraction. The results are listed in Table 01 We note the mean values, n^, for the charged 
and neutral modes are consistent within errors to zero signal events. 

The 90% confidence level (C.L.) upper limit number of events (A'^go) is calculated using the 
Gaussian probability distribution with the assumption the number of signal events is always greater 
than zero. The integral of the distribution from zero to Nqq will be 90% of the total area above 
zero. Combining Nqq, e, N^g events, and the secondary branching fractions, we obtain. 
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5.8 X 10-*^ (90% C.L.) 



eNjj^B{J/i) l+l-)B{K^s vr+vr-) 
11 X 10~^ (90% C.L.). 



for the neutral and charged branching fraction upper limits. For completeness we include the 
central value (68% confidence interval) for the branching fraction using the Ug ± Gs values, 



(-5.5 ±5.2) X 10-^ 



(2.3 ±5.7) X 10"^ 



The results are summarized in Table 13 

Table 3: The estimated number of signal events, 90% C.L. upper limit of signal events, the branching 
fraction upper limits, and the branching fraction (B) for each i?-decay. 



Mode Us ± as Ngp 90% C.L. B 

J/iIj-k'^tt^K^ -14.8 ±13.9 15.6 < 5.8 x lO"'^ (-5.5 ± 5.2) x 10"' 
J/i>-R^Tv°K'^ 3.4 ±8.3 15.9 < 11 x 10"^ (2.3 ± 5.7) x 10"^ 



5 PHYSICS INTERPRETATION 

We test the charged partner hypothesis at a mass of 3872 MeV/c^ using a likelihood ratio 
test |15j . Here we determine the ratio of the two probabilities from the null (Hq) and signal {Hi) 
hypotheses using our experimental observation of 87 events in the signal-box. 

The null hypothesis assumes the estimated background events, rihzizab, produced all the observed 
signal-box events. Assuming the background probability distribution is a Gaussian function, we 
calculate a probability of P{Hq)=7.34: x 10^^ to measure 87 or fewer events. 

The isovector signal hypothesis predicts the product branching fractions to have the ratio, 
B{B X^K,X^ J/i^P^) =2 B{B X(3872)K, X(3872) ^ J/V'P°)- Using the BABAR 
branching fractions B{B^ X{5872)K^, X(3872) J/V'tt+tt^) = (1.28 ± .41) x 10"^ and 
assuming all the vr+vr^ decays originate from p*^, we expect B{B X^K^, X^ J/tpp^) = 
(2.56 ± 0.82) X 10^^. This would produce 69 ± 23 observed signal events in a data sample of 
213 million BB events. The error combines the uncertainty on the branching fraction and the 
systematic error, asys, on our efficiency. The probability distributions for the signal events and the 
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estimated background events are modeled as two uncorrelated Gaussian functions. The probability 
of observing 87 or fewer events with this probability distribution is P(-ffi)=1.18 x 10~^. 

The likelihood ratio (A) test of the null hypothesis relative to the signal hypothesis yields A = 
P{Hq) / P{Hi) = 622. This corresponds to a probability of less than 1 part in 600 that the 
hypothesis is correct with the outcome of our measurement. Hence our result does not support the 
existence of charged molecular states or charged partners of the X(3872). 



6 SUMMARY 

In conclusion, we have performed a search for a charged partner of the X(3872) decaying to 
J/i/jn^TT^. Our results set upper limits on the product branching fractions of B{B^ / X^K^ , 
X^ J/il^n^-K^) < 5.2 X 10-^ and B{B^ X^K^, ^ J/ip'K^n^) < 11 x 10"^ at the 90% 
confidence level. We exclude the isovector X hypothesis with a likelihood ratio test and with 
our experimental results we obtain a ratio greater than 600 for the null hypothesis relative to the 
isovector signal hypothesis. 
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Figure 1: The /S.E (a) and mss (b) distributions from the B^/B^ — > J/i/jit^tt^K^ mode after 
applying the optimized cuts. 




Figure 2: The AE (a) and mES (b) distributions from the 5=*= — J/ip-ir^Tr'^Kg mode after applying 
the optimized cuts. 
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Figure 3: The tt^tt^K'^ (a) distributions from the ^ J/V'7r=^7r°K=F mode and the tt^tt'^K^ 

(b) distributions from the J/iptr^Tr^Kg mode after a signal-box and a mass cut. The 

dashed histogram is obtained from events in a 5.24< ttt-es <5.26 GeV/c^ sideband. The neutral 
and charged i^i(1270) signals are evident in the plots (a) and (b), respectively. 




m(J/'J'it-ic'') (GeV/c^) m(J/'J'it-7c'') (GeV/c^) 

Figure 4: The J/ipTr^iT^ invariant mass from neutral (a), B^/B^ —>■ J /ijjTT^iT^K^ , and charged (b), 
B^ J/tPTT^TT^K'^, modes. 
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Figure 5: The J/ip-K^Ti^ invariant mass in 10 MeV/c^ bins for the neutral (a), B^/B^ — > 
J/il^TT^TT^K^, and charged (b), B^ J/iI^tt^tt^K^, modes. These plots are the same as the 
previous Figs, except they have finer binning. No evidence for the charged X(3872) partner is 
observed in either plot. 




Figure 6: Fitted tties distribution for the B^ mode with the signal region selection, 
- 3872MeV/c2| < 12MeV/c2, (a) and with the sideband region selection, 48 < 
|m( J/V^vr^vr'^) — 3872MeV/c^| < 72MeV/c^, (b). The signal region selection plot is used to es- 
timate the combinatoric background. The sideband region selection plot is used to estimate the 
peaking background. The sideband region selection has twice the J/ipTr^ir^ mass range of the signal 
region selection. 
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Figure 7: Fitted rriES distribution for the mode with the signal selection, \m{J/ipiT^7r^) — 
3872MeV/c2| < l2MeV/c^, (a) and with the sideband region selection, 48 < |m(J/'i/'vr^7r'^) — 
3872 MeV/ < 72 MeV/ , (b) . The signal region selection plot is used to estimate the combinatoric 
background. The sideband region selection plot is used to estimate the peaking background. The 
sideband region selection has twice the J/ipTT^ir^ mass range of the signal region selection. 
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